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Abstract 

We study a simplified model of the RNA molecule proposed by G. Vernizzi, 
H. Orland and A. Zee in the regime of strong concentration of positive ions in 
solution. The model considers a flexible chain of equal bases that can pairwise 
interact with any other one along the chain, while preserving the property 
of saturation of the interactions. In the regime considered, we observe the 
emergence of a critical temperature T c separating two phases that can be 
characterized by the topology of the predominant configurations: in the large 
temperature regime, the dominant configurations of the molecule have very 
large genera (of the order of the size of the molecule), corresponding to a 
complex topology, whereas in the opposite regime of low temperatures, the 
dominant configurations are simple and have the topology of a sphere. We 
determine that this topological phase transition is of first order and provide 
an analytic expression for T c . The regime studied for this model exhibits 
analogies with that for the dense polymer systems studied by de Gennes. 

PACS numbers: 87.14gn , 02.10.Yn , 87.15. Cc 
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1 Introduction 



Recent advances in mathematical models of RNA and DNA molecules have led to 
stimulating new studies of their complex spatial properties and their relation with 
the environmental variables ( see, e.g. [H El El HI E] ). In this paper we consider 
the exactly solved model for a RNA molecule presented in [6], which consists of a 
homopolymer of length L, with an infinitely flexible backbone, and in which any 
arbitrary pair of bases is allowed to (pairwise) interact. This combinatorial model 
preserves the important property of saturation of the interactions [21 [7] of the actual 
RNA molecule, but does not include both the geometric and energetic aspects of 
the real molecules. Due to its simplicity, the model is exactly solved, i.e., it allows 
for analytical expressions of statistical and topological interesting quantities, such 
as the partition function and its topological expansion [§]. A crucial ingredient of 
the model is the introduction of an extra degree of freedom, N, such that a random 
N x N hermitian matrix is added to each base position along the chain. A physical 
interpretation can be given to this parameter a posteriori 0,181 [9]: 1/iV plays a role in 
the model that can be associated to the concentration of positive ions (like Mg ++ ) 
in solution for real molecules. These positive ions are responsible for the overall 
electric charge neutralization of the system compensating the negative charge in 
the phosphate groups in the RNA molecule. The neutralization of the phosphate 
ions provide the physical mechanism for the folding of the molecule [TOl [TT1 [T2l 
H31 HU 02] ■ Moreover, the algebraic power of 1/N 2 in the topological expansion of 
the partition function of the model is the genus of the diagrams associated to the 
different configurations [El El HE]. These diagrams encode the real spatial structure 
of the molecule, in the sense that they are analogous to the secondary (or planar) 
structure of the molecule, from which the tertiary (or spatial) one is obtained [6j. 
Therefore, the parameter N works gulator that controls the topology of the 

molecule within the model. To the best of our knowledge, the regime of N — > 
has not been received much attention in the theoretical models considered in the 
literature, and in this paper we address some questions arising in this regime. We 
identify this regime with the dense polymer phase of the 0(n) vector model of de 
Gennes [17J, which is described by the limit n — > 0. Furthermore, we are not aware 
of any experimental works in real RNA molecules in the corresponding laboratory 
regime. 

In a previous work |B], we have provided an exact analytic expression for the 
partition function of the model, which has been used to study its thermodynamic 
properties in the regime N > 1, i.e., that of small ion concentration in the sur- 
rounding environment. However, some of the expressions that we have obtained in 
that paper are analytical in N and are, therefore, useful to explore other regimes 
than the one we have chosen to discuss in that work. In this paper we study (for a 
wide range of temperatures) the complementary regime N — > , which can be char- 
acterized as one with large positive ion concentration in the surrounding medium. 
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We observe the emergence of a critical temperature T c separating two regimes with 
different spatial configurations: in the large temperature regime, the dominant con- 
figurations of the molecule have very large genera (of the order of the size of the 
molecule), corresponding to a complex topology, whereas in the opposite regime of 
low temperatures, the dominant configurations are simple and have the topology of a 
sphere. This transition is not related a-priori to the so-called coil-globule transition 
in which the molecule adopts the shape of an elongated coil for low temperatures, 
whereas it assumes a compact globule shape in the large temperature regime. This 
latter transition is considered to be of first order 0, HH1 EES EH EI] . 



Some recent work in the literature has recently addressed related issues. In [18J, 

a study of the distribution of genera of pseudoknotted configurations in the In- 
dependent phase transition for a self-avoiding homopolymer on a lattice has been 
considered. A study of the behavior of the distribution of genera for the model 
introduced in p] with external perturbations has been the subject of Refs. [221 EH]- 
Moreover, in [22], a 'structural transition' dependent on the external perturbation 
has been also discussed. 

This paper is organized as follows: in Section 2 we reconsider the model pre- 
sented in [U] and recall some exact results, such as the partition function and its 
topological expansion [8]. Section 3 is devoted to the study of the regime N — > 
and the topological phase transition: we first consider the mathematical setting and 
approximations that are needed to analytically treat the phase transition. We then 
apply this preliminary study to the partition function of the model to make explicit 
the emergence of the phase transition and we give an analytic expression for the 
critical temperature. In section 4, we discuss the thermodynamic properties of the 
system in the limit N —>■ 0. We calculate the free energy per particle, the entropy 
and the latent heat. From the last thermodynamic quantity, we show that the phase 
transition is of first order. Finally, we give our conclusions. 

2 Topological expansion of the partition function 

We start by reviewing some of the properties of the model proposed by G. Vernizzi, 
H. Orland and A. Zee [6]. The model considers a chain of bases of one type only, 
such that the interaction energy between any pair of bases is a constant e. A given 
base can interact with any other base in the chain, but preserving the saturation 
property, which excludes interactions among three or more bases. We consider that 
this property is one of the essential ingredients of the simplified model. Therefore, 
all the Boltzmann factors v = exp (— e//«T) (where T is the absolute temperature 
and k is the Boltzmann constant) are equal as well. At each base site, a random 
N x N hermitian matrix is added as the relevant degree of freedom. We consider 
this feature to be a second essential ingredient of the model. The configurational 
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partition function Z of a chain of length L is: 

Z(L,N,T) = I(tr(l + ^) L ) (2.1) 

d<p e-^^ 2 ^tr(l + f dtp e~^\ (2.2) 

where (p represents a collective degree of freedom (a sort of center-of-mass random 
matrix) . The simple form of (12.21) is a consequence of the symmetry of the matrix 
potential that reduces the original integration over L matrices to one integration 
over p [6J. Applying standard results of random matrices to f)2.2p one obtains: 

[L/2] 

Z(L, N,T) = J2 ML, N) e-^ T , (2.3) 

fc=0 

where the symbol [L/2] means the integer part of L/2, = ke and 

«^-£U)(i)(A)*^- (2 - 4) 

From (12. 3p we may compute Z exactly (for each L), as a function of A" and T. The 
spectrum of the system has [L/2] + 1 energy levels, with energies: 0, e, 2e, . . . , [L/2]e 
and the degeneracy of the k — th level is dk(L, 1). 

As we have mentioned before, the power of 1/A^ 2 yields the genus g of the diagram 
obtained from the chain by joining all interacting bases by a 'photonic line' [BJ, that 
is, the minimum number of handles of the surface on which the diagram can be 
drawn without crossings. In [S], we have written the partition function of the model 
in the form of a topological expansion [TBJ El E],ie., as a power series in l/N 2 , 
where the coefficients take into account all the Feynman diagrams with the same 
topological character: 

Z(L,N,T) = J2z 9 (L>T)— . (2.5) 

9=0 

Here z g (L, oo) is the number of planar diagrams that can be drawn on a topological 
surface of genus g for a molecule of size L. Note that, as a function of T, z g (L,T) 
is the partition function of the system living on the topological surface of genus g. 
The coefficients z g (L,T) are given by: 

[L/2] k 2j-k g(fc+l-2<?) 

where is the Stirling number of the first kind [2U [25] with parameters m,j 

(5'- m ^ = if m > j or if j < 0). In the limit T — > oo, z g (L,T) coincides with the 
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coefficients a>L,g of Ref. [6j. Using the property of the Stirling numbers mentioned 
in this paragraph, we see from (I2.6P that the maximum genus of a diagram for a 
given L is [L/4] and, therefore, g < [L/4]. 



3 The partition function in the limit N —> 

3.1 Preliminary study of the analytic behavior of the par- 
tition function 

In a previous paper, we have considered the limit N — > oo of the partition function 
( 12.31) [H] , for which we gave explicit expressions. As a consequence, we have verified 
the consistency of interpreting the parameter 1/JV as being proportional to the den- 
sity of positive ions in the media surrounding the molecule, as has been suggested in 
[U] . Therefore, the range of the analysis considered in [S] applies to media with small 
concentration of positive ions. It is natural to consider as well the opposite situation 
of large positive ion concentration in the surrounding medium, which corresponds 
to the limit iV — > 0. We are going to show later on that this regime is characterized 
by the existence of a critical temperature and a phase transition. 

From the mathematical point of view, studying this limit requires a careful han- 
dling of the partition function considered as an analytic function of several variables. 
Here we discuss a mathematical scheme that allows us to obtain simple analytic ex- 
pressions in this limit. We first consider a general function of two variables N and 
T (which will be ultimately related to the free energy) which admits a decomposition 
of the form: 



where 0/ (/ means 'fast') indicates a function whose asymptotic grow in the vari- 
able T is in the class of the exponential function, whereas the function S (s 
means 'slow') denotes a function that grows slower than the class of the exponen- 
tial function (e.g., polynomial, logarithmic classes). Depending on the value of N 
(which we consider fixed at Nq) there exists a small range of values of T for which 
0(<f)f(T)) = 0((f> s (No,T)) for values of T nearby a critical value T c . At this point 
we can anticipate the idea of defining the critical parameter T c as the value of T for 
which both the slow and fast parts of the function become equal. We shall further 
develop on this idea below. The parameter T c emerges naturally in the case when 
the function is the free energy where it plays the role of a critical temperature. 
Given the asymptotic behavior of 0/, as soon as T crosses T c , the order of 0/ varies 
rapidly, increasing or decreasing, depending of the sense in which T crosses T c , so 
that (for example): 



0(iV,T)=0 / (T) + 0,(iV,T) , 



(3.1) 




(3.2) 
(3.3) 
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Therefore, <f)(N , T) has different analytical behavior according to the region in which 
T is located with respect to T c . The change in between these two regions becomes 
more pronounced when <pf varies more rapidly. In Fig. [T]we show an example of 
this behavior for the cases 0/ = e _1//T and <p s = NT 5 (in this example T and N 
are two dimensionless variables) , in the vicinity of T c ~ 0.00466 (the plot has been 
done numerically): for T > T c , is dominated by 0/ and for T < T c it is dominated 
by 0s- 
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3. x 1(T 92 
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0.0045 ~ T c 


0.0050 



Figure 1: Plot of 0(iV, T) = e' 1 ^ +NT 5 around of T c ~ 0.00466 with N = 10- 80 .The 
plot has been done numerically. 

If we assume that <pf is of the form 0/(T) ~ e _1//T , the change is even more 
noticeable for small values of T . If we had T c — > 0, the derivative <J)'(Nq,T = T c ) 
would not exist. We shall correlate the absence of the derivative of at the critical 
point with the appearance of a phase transition of the system for T c = 0. 



3.2 The critical temperature and the phase transition 

We now consider the topological expansion in 1/N 2 (12.51) of the partition function Z 
of the system. Considering that the maximum genus of the diagrams of all possible 
configurations is [L/4] (as we have mentioned above), we rewrite Z as a polynomial 
in the variable iV 2 times a iV-dependent factor which is divergent in the N — ► 
limit: 

Z(L,N,T) = ^^Z fl (L,iV,T), (3.4) 
Z R (L,N,T) = Z[m + Z[L/4]-iiV 2 + ■ ■ ■ + VV 2[L/4] . (3.5) 
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Here the subindex R indicates that the function is regular in the limit N — > 0. The 
coefficients z g were given in the previous section [8], and can be written as: 



[L/2] 



2 </ 



(^^) = E r ^( L ) e " fce/KT ' ( 3 - 6 ) 



fc=0 

where the coefficients rf. :9 {L) are given by 

r ^ (L) = (L-2fc)!2* . = 5 2fl •./)!/!!./ • 1)! ' (3 ' 7) 

and Sl j) are the Stirling numbers of the first kind with parameters i and j. Using 
the property = 5io [26], we rewrite the partition function on the topological 
surface of genus [L/4] — n as: 

[L/2] 

z im -n(L,T) = Yl r k,[L/4]-n{L) e~ he/KT , (3.8) 

fc=[L/2]-2n+<5 

where 5 = — 1 for [L/2] odd and zero otherwise. Without loss of generality, we shall 
consider [L/2] even from now on. Substituting the coefficients z g in ( 13. 51) . we obtain 
the regular part of the partition function (RZ): 

Z R (L, N, T) = e-^ KT {r [L/2Um + {r [L/2] ^ [m ^e 2 ^ T + r [m _ um _^ KT 

+ r [i/ 2], [ L/4]-i)iV 2 + ... + (r , e^ 2 ]^ 

+ rii0e ([V2]-i),/«T + . . . + r[L/2])0 )^[x/4] } (3 g) 

Note that r^o = 1 since it counts the number of configurations with lowest energy 
and genus, that is, the unique configuration without interaction among bases. In 
order to treat (13. 9p analytically in what follows, we shall consider two approximations 
on it, valid in the limit N <C 1. We have studied the analytic behavior of (13. 9p 
using the program Mathematica [27J, within the range N J$ 10~ 20 in the concrete 
examples considered. For these examples, the range of sizes of the system is L ~ 
10 2 — 10 4 and that for the temperature is T ~ 10~ 4 — 10~ 2 . 

First, we observe that for the ranges of L and T considered, the terms proportional 
to N 2g are dominated by exponential factors {e.g., r k)g /r k _i g ;$ 10 7 and e^r > 10 43 
(e = k — 1)). For large values of T the most important term in the sum (13.91) is the 
first one (up to order zero in N), whereas for small values of T, the dominant term 
in the sum is the one proportional to jV 2 ' L//4 l (= jV' L//2 l, for [L/2] even). To keep 
analytic expressions as simple as possible, we retain only the most important term 
proportional to N^ L ^ 2 ^ in (|3.9|) . i.e., e^l 2 ^! 1 ^ and write: 

Z R {L, N, T) ~ r [Lmm e-^'* T + . (3.10) 
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We refer to f 1 3 . X j) as the Thermal Approximation (TA) for the partition function, 
since the relevance of the terms in Zr vary significatively when T crosses T c , that 
is, when T satisfies 0(r [L/2UL/4] e^ L / 2 ^ KT ) = O(N^) (see EI]). 

Next, we specify the definition of T c with more detail. It is clear that (13.101) is 
a broader approximation to (13.91) . since it excludes many terms that are relatively 
important when the two terms in (13. lOf) are of the same order, i.e., when T is closer 
to T c . However, we will see later how to minimize the difference between ( 13. ldj) and 
(13.91) . showing the conditions under which RZ becomes close to TA. From the TA 
we calculate the regular part of free energy from F = — kT In Z: 

F R (L,N,T) ~ - K Tln{r [L/2]i[L/4] e-[ L / 2 ]^ T + iV[ L / 2 ]} . (3.11) 

In Fig. [5] we show a plot of the thermal approximation for the free energy TA against 
the temperature. Two regimes can be clearly differentiated, separated by T = T c : 
for T < T c , the second term in the logarithm function of (13.111) is much larger than 
the first one, and can be therefore neglected, with opposite behavior for T > T c . 
These plots show two interesting characteristics: (13. lOj) behaves linearly below and 
above T c and the slope for T < T c is higher than that for the case with T > T c . We 
therefore identify a phase transition occurring for T = T c , which we identify as the 
critical temperature. As it is customary for a finite system such as the model for 
the RNA molecule studied here, the term 'phase transition' should be interpreted 
'strongly cooperative phenomenon' |28j . 




Figure 2: Typical behavior of Fr calculated from RZ and TA as a function of T for 
L = 100 and N = 10~ 120 . 

In a real RNA molecule, the large negative electric charge of the phosphate ions, 
prevents it to fold onto a compact structure. The concentration of positive ions in 
solution, such as Mg ++ , neutralize the phosphate ions making possible the folding 
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of the RNA molecule [291 DSl EE]- In the context of the model studied in this paper, 
if we interpret 1/N as the concentration of Mg ++ in solution [El [21 E], we see from 
(13. lip that increasing this concentration favors the formation of RNA structures 
with large genus, and viceversa. 

The asymptotic straight lines in Fig. [2] are obtained by retaining only the leading 
term in the logarithm of (I3.1ip when varying the temperature around T c . Therefore, 
the analytic expressions for these asymptotic straight lines are given by: 

Fr (L, N,T < T c ) ~ -K[L/2]ln{N}T , (3.12) 
F R {L ) N 1 T>T C ) ~ -Kln{r [i/2]i[L/4] }T+[L/2]e, (3.13) 

and we obtain an analytic expression for the critical temperature equating the two 
terms in the logarithm of (13.111) : 

T C (L, AO = - 1 — T (3.14) 

Kin \(r [L/2Um )^jj 



As we have mentioned in Sec. 13.11 and show in Fig. [21 the regular free energy 
Fr undergoes a change in its analytic behavior when the temperature reaches the 
critical value T c . Below T c , Fr is dominated by the slowly varying function (N^), 

whereas above T c , the rapidly varying function (rr£/2],[z/4]e~' L «T~) dominates. The 
change in the behavior of Fr is sharp in the TA, while for RZ it is obtained for 
a small range of temperatures because RZ presents additional terms that acquire 
relative importance when T is close to T c . 

It can be seen numerically from (13. 7j) that the function of L in T c behaves linearly 

i 

as {t[l/2}\l/a}) [l/2] = L/e — 4.5521. For large L (L > 200), we can write the critical 
temperature as: 

(315) 

Note that (13.151) possess the symmetry T C (L, N) = T c (aL, aN) for a > a real 
parameter, that is, T c is scale invariant. Therefore, in the large L limit, we have 
that T C (L, N) = T C (L/N) Furthermore, there exists a natural cut-off for iV such that 
T c > 0: 

N < (r[i/2],[L/4])W ~ L/e . (3.16) 

This is a consequence of (13. 14ft and it is consistent with the condition stating that 
the slope of the asymptotic straight line for T < T c must be larger than that for 
T > T c . However, this cut-off is not very restrictive because we are interested in the 
the limit iV 1 and L is at least of 0(1) . 

Moreover, there exist further and more restrictive conditions on N, implied by 
both the thermodynamic and large size limits. Given that RZ should approach TA 
for large L, a relationship between iV and L should exist. When L increases and N is 
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fixed, the difference between the analytic expressions for RZ and TA also increases. 
In order to obtain a convergence between the two in the thermodynamic limit, we 
make iV a L-dependent variable as follows: for each L, there exists an upper limit 
for iV below which the phase transition exists. This can be seen by demanding that 
the scenario discussed in the previous section actually occurs. For real molecules, 
this would mean that there should exist a minimal (L-dependent) concentration of 
Mg ++ for which the phase transition exists. As a consequence of this restriction, in 
the thermodynamic limit L — > oo then N — > , implying T c — > 0. In this limit, the 
derivative of Fr at T — T c (see (13. IIP ) does not exist and there is a bona fide phase 
transition. 

Furthermore, the limit N — > in the model we have considered can be associated 
with the dense polymer phase of the n vector model which arises for n — > as has 
been discussed by de Gennes [17J. The correspondence is established between the 
degrees of freedom of the 0(n) vector model, which are n-component spin vectors, 
and the degrees of freedom of the RNA model of [6j, given by by N x iV hermitian 
matrices. The limit n — > in the de Gennes 0(n) model corresponds to a high 
density polymer phase, which is naturally associated to the high concentration of 
Mg ++ (or other positive ions such as K + ) in solution phase of the RNA model. 

From the previous discussion of the RNA model, the T-dependent phase transition 
involves a topology change in the spatial configurations of the molecule, which goes 
from one with genus zero for T < T c to another one with large genus [L/4] for 
T > T c . The difference in the topology of the configurations of the molecule when T 
crosses T c justifies the use of the term 'topological' for describing the nature of the 
phase transition. Note that this transition does not correspond to the coil-globule 
transition in polymers studied in [T5], in which the genus of the configurations 
decreases with increasing temperature (see Fig. [3]). 



high[Mg ++ ] 
low [Mg ++ ] 

Globule state 




Coil state 



Figure 3: Phase transition between coil and globule states when T crosses T c . The 
arrows indicate the direction in which T increases, to right (left) for high (low) 
concentration of Mg ++ . This is not to confused with the topological phase transition, 
as explained in the text. 

This is due to the fact that both transitions take place in different regimes: 
whereas in [18J the RNA molecule is surrounded by a dilute solution of Mg ++ , 
for the case we have discussed the molecule is imbedded in a medium with large 
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concentration of positive ions. 



4 Thermodynamic properties in the N —> limit 

From the analytic results for the partition function obtained in the previous section, 
we now calculate some thermodynamic quantities that we will use to better describe 
the configurations of the RNA molecule above and bellow of T c and characterize the 
nature of the transition that it undergoes. It is important to remark here that the 
condition iV < 1 (interpreted as meaning a large concentration of Mg ++ ) is valid 
for all the results presented in this section, since this is a necessary condition for 
using the approximate expression for the partition function instead of the exact one 
(13.91) . In this analysis, we only consider the regular part of the partition function, 
disregarding the iV-dependent, divergent factor in (13.51) . This factor can be absorbed 
by convenient renormalizations, and cancels out in any observable quantity which 
involves an statistical average over the ensemble of random configurations of the 
system. Therefore, all thermodynamic quantities studied will be, therefore, termed 
as 'regular'. 

4.1 Free energy per base 

From (13.101) . we determine the free energy per base, / = F/L. Using that [L/2]/L — > 
1/2 in the large L limit, we write the regular free energy per base as: 

f R (L,N,T<T c ) ~ -~HN}T (4.1) 

f R (L,N,T>T c ) ~ _^ln{L/e}T + | . (4.2) 

On the one hand, from (13 . 12[) we notice that /# is independent of L for T < T c . 
This means that f R is determined by the local (short distances) environment of a 
point in the chain and not by global (large distances) properties. This behavior is 
reasonable because thermal agitation is small, which prevents coupling on a given 
base with another non-neighboring one in the chain. On the other hand, for T > T c 
we see from ( 13 . 1 3f) that fn is in turn independent of N. This can be interpreted 
by noting that for temperatures above T c , thermal agitation overcomes the folding 
action of the positive ion concentration in the medium, rendering the local behavior 
independent of this concentration, and, therefore, of N. 

4.2 Entropy 

Next, we calculate the regular entropy Sr = —OFr/ST from RZ and TA, and plot 
it in Fig. H] , which shows that the dependence of the entropy with the temperature 
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resembles a step function, with a positive step in T = T c (in particular, for the TA 
case). This behavior could be expected, since Sji(T a ) < Sjj(T 6 ) for T a < T b and it is 
a consequence of the linear asymptotic behavior of the regular free energy displayed 
in Fig. [2j On the contrary, a negative step in T = T c would have implied T c < 0. 
The curves of Sr from the TA, for T below and above T c , are derived from f l 3 . X 2 f) and 
(I3.13p . Moreover, it can be seen from Fig. HJthat as N decreases, RZ approaches 
TA and the behavior of RZ resembles more closely that of a step function: 

S R (L,N,T<T C ) ~ 4L/2}\n{N} , (4.3) 
S R (L,N,T>T C ) ~ K\n{r [L/2Um }~K[L/2]\n{L/e} . (4.4) 

We have verified numerically that, in the limit iV — - > 0, the plot of Sr from RZ 
behaves asymptotically as step function with step in T c of magnitude proportional 
to the latent heat. 




Figure 4: Typical behavior of Sr vs. T calculated from RZ (continuous line) and 
TA (dashed line), for L = 200, temperatures close to T c and different values of iV: 
(a) N = 10~ 80 , (b) N = lO" 100 , (c) N = 10~ 130 , (d) = lO" 180 . The convergence 
of RZ to TA when N decreases is displayed. 
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4.3 Latent heat 



From the expressions for the entropy, we calculate the regular latent heat of the 
transition: 

l R (L,N,T = T c ) = T C AS R (T C ) = T C {S R (T > T c ) - S R (T < T c )} = [L/2\e . (4.5) 

Eq. (14. 5 p shows that the phase transition in the model we have considered is of first 
order [3U]. Note that (14. 5p does not depend on N: this result could be expected, 
since the latent heat is the energy released or absorbed by the system during the 
transition, and N is associated with the concentration of Mg ++ , and, therefore, plays 
the effective role of an external variable regardless of the way in which it has been 
introduced. Eq. ( 14.51) expresses that the energy exchange between the system and 
the bath during the phase transition is [L/2]e, since there are [L/2] parings between 
pairs of bases, which could be created or broken during the transition, depending 
on the direction in which T crosses T c . 

5 Conclusions 

In this paper we have studied, using both analytical and numerical methods, the 
T— dependent phase transition in a simplified model of the RNA molecule, in the 
regime of large concentration of positive ions in solution. This regime has similarities 
with the large density phase of polymers studied by de Gennes. We have presented 
an analytical expression for the critical temperature T c , which tends to zero in the 
thermodynamic limit. The critical temperature separates the only configuration 
without interaction between the bases (and, therefore, of genus zero) which is dom- 
inant for T < T c , from the configurations with high energy and large genera equal 
to [L/4], which dominate for T > T c . This transition is no to be confused with 
the coil-globule transition, which appears for low concentration of positive ions in 
solution. Due to the interesting dependence of the genus with the temperature, we 
call this a topological phase transition, and we have shown that the transition is of 
first order. Although the model studied is very simple, we are confident that most 
of the properties studied here might be robust and extend to more realistic ones, 
given the importance of the saturation property of the elementary base interaction, 
that this model preserves. 
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